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Abstract 

We consider the exponential moments of integrated currents of ID asymmetric 
simple exclusion process using the duality found by Schiitz. For the ASEP on the 
infinite lattice we show that the nth moment is reduced to the problem of the ASEP 
with less than or equal to n particles. 

1 Introduction 

The one dimensional asymmetric simple exclusion process (ASEP) is a many-particle 
stochastic process in which each particle is an asymmetric random walker but with ex- 
clusion interaction among particles [9j|T0] . We consider the ASEP with hopping rate p to 
the right and q to the left with p + q = 1. The ratio is denoted by r = p/q. 

The ASEP can be denned either on a finite lattice or on an infinite lattice. It is of 
much current interest to study fluctuation properties of the ASEP on Z because for this 
case one can perform detailed analysis using the connection to random matrix theory and 
other techniques [31 [HI EH E2]- Until a few years ago the analysis had been restricted to 
the totally asymmetric case, p = or q = 0, i.e., particles hop only in one direction. 

In [25], Tracy and Widom succeeded in computing the distribution of the particle po- 
sition for the ASEP with general parameter values using the transition probability derived 
from the Bethe ansatz. For recent developments see P,E31EHEErEI] ■ 

Recently their formula was utilized to study the height fluctuations of the KPZ equation 
[^U5yT5HTT] . By the Cole-Hopf transformation, the KPZ equation is mapped to a problem of 
directed polymer. As noted long time ago, the n replica partition function of this directed 
polymer is mapped to the problem of attractive 5-Bose gas with n particles [Hill]. This 
connection has been utilized recently in [TjE] . By considering a generating function of the 
n replicas, one could reproduce the reslults for height fluctuations of the KPZ eqaution. 
See also a related work \12\. 
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In this note we point out that similar consideration is possible for the ASEP. We show 
that the nth exponential moments of current in the ASEP can be written as a summation 
of transition probabilities for k particles with < k < n of the ASEP. From this one could 
find an expression for current fluctuations of the ASEP. Our argument is based on the 
duality relation of the ASEP found by Schiitz [18] . 



2 Duality 

First we consider the ASEP on a finite lattice [L] = {1,2, ... ,L} with reflective boundaries. 
We employ the formulation using the quantum spin chain language, see e.g. |20j. We set 
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(2.1) 



Let us introduce a vector |0) which corresponds to the empty system. One can construct 
an n particle state with particle positions at x\, . . . , x n by 

\x 1 ,...,x n ) = s~ l ---s-JO). (2.2) 

Here s~ means it acts nontrivially only on the space of site x as a 2 x 2 matrix s~ in (12. ip . 
The state of the system can also be specified by a set of particle numbers rj x on each site 
x G [L]. Here rj x — 1 (resp. rj x = 0) means the site x is occupied (resp. empty). We 
sometimes abbreviate rj = {r/i, . . . ,t]l} and denote the corresponding state by \r)). 

Let P(r], t) be the probability that the configuration of the system is rj at time t and 

set 

\P) = J2P(v,t)\v) (2.3) 
v 

where Yl v means the summation over all particle configuration. The time evolution of this 
is given by the master equation, 



d_ 

di 



P) = -H\P), 



(2.4) 



where the transition rate matrix is given by 

L-l 

H = ~ [P s t s J+i + <? s 7 s i+i - Pn>j(l - n j+1 ) - q(l - nj)n j+l ] . 
i=i 

The transition probability, i.e., the probability that n particles starting from yi, 
time are on sites at time t is written as 



G(xi, ...,x n ; t\yi, . . . , y n ; 0) = (xi, . . . , x n \e 



-tHi 



yi, 



Vn) 



(2.5) 
y n at 

(2.6) 



with (xi, . . . , x r , 



the understanding {x} n 



s^ n . We sometimes abbreviate this as G({x} n ; t\{y} n ; 0) with 
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We recall the duality of the ASEP based on a quantum group symmetry of the process 
. Our notation in this article is slightly different from [TS]. Let us set 



L 



X + = T ~L/4 + l ^ r E*= 1 1 d-n J ) s - j (2 .7) 



fc=l 

L 



k=l 

# = T -E*=i»J. (2.9) 
They satisfy the U q (sl2) algebra 



KX+R- 1 = tX + , KX~K~ l = r^X-, (2.10) 

K-K- 

rV2 _ T - 



^ + ^1 = ,i/ 2 _ r -i/2 - ( 2 - n ^ 



They commute with H, i.e., 

[H,X ± } = [H,K]=0. (2.12) 

The last relation reflects the fact that the number of particles is conserved in the dynamics. 
Since the state space is finite there is the unique stationary state for each particle number 
N,0 < N < L. Let us denote it by \N). It satisfies H\N) = 0. When N = this is 
nothing but the state |0) with no particle. For N > 1, one can construct \N) by applying 
X + for iV times to |0) as 

\N) = (X + f\0). (2.13) 
It is easy to check H\N) = using (12. 12)) and H\0) = 0. If we introduce 

(N\=C N (0\(X-f = Yl W ( 2 - 14 ) 



where 



LN 



Cn = t~- ^ 7^ T/vT' (2-15) 



'1-T 1-T 



the normalization of a state \P) is written as (N\P) = 1. Using (12.71) . one sees that the 
normalized version of (12.131) is given by 

m^ = r-^ ^^-U S ^ ,t "* l " > - (2 - 16) 

V-I2k=iVk=N 

The average of a quantity A which depends on rji's at time £ for an initial state |ijy) with 
N particle is 

(A) t = (N\Ae- Ht \I N ). (2.17) 
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Now let us define 



and set 



iV. = X>i (2-18) 



One can verify 
Lemma 1. 

where 



Qx = t"% (2.19) 
4 = Qx Z ^ = t*- 1 ^. (2-20) 

T — 1 



l(X-) N , Q x ] = (r N - l)Q x X-{X-) N - 1 (2.21) 



X- = r-^^s+r^+i™,. (2.22) 
fc=i 

Proof of Lemma [ZJ The N — 1 case is easily checked by using the following relations 

T n s + = s+, s~r n = rs + , r n S - = ts~, s~r n = s~ . (2.23) 

Next assume (I2.2ip is true for N. We want to see (12.211) holds for AT + 1. We start from 
[{X-) N+1 ,Q m ] = l(X-) N ,Q x ]X- + (X-) N [X-,Q x }. (2.24) 

One uses (T2T2TD for A?" = 1 and N to get 

[{x-) N+ \Q*) 

= (r N - l)Q x X-(X- f + (t n - 1)( T - l)Q x X-{X-) N - x X- + ( T - l)g,pr)"Jf- 

(2.25) 

Comparing this with RHS of f)2.2ip for AT + 1 , it is enough to show 

(X- - X-)(X-f- l X- = r N X-(X-) N (X- - X-). (2.26) 

To verify this, one can check the AT = 1 case by using (I2.8p . fl2.23p and then use mathemat- 
ical induction. 

□ 

Applying this lemma, we see, when distinct, 

(N\Q X1 ■■■Q Xn = C N , n ( Xl , . . . , x n \(X-f- n (2.27) 

where 

T \(N-n)LM _ T ^N~n 

° N > n = (l-r)-(l-r^)- (2 ' 28) 
Notice Cjv = Cjv,o- Using this one can show 



Proposition 2 ( [IB])- When Xi's are distinct, it holds 

{Qx 1 ---Qxjt= Yl G (xi,---, x n;t\y 1 ,...,y n ;0)(Q yi ---Q yn )o- (2.29) 

1<2/1<— <Vn<L 

Proof of Proposition This is seen as follows: 

mS=(N\Q xl ---Q Xn e- Ht \I N ) 

\(X-) N - n e- Ht \I N ) 

.71 \ X 1 5 • • • 5 X 7l 

\e- H \X-) N - n \I N ) 
= ^ ( x i'-- -i x n\e~~ Ht \yi, ■■■iVn)- C N)n {y u ... ,y n \(X-) N ~ n \I N ) 

1<2/1<— <Vn<L 

(xi,-.., x n \e- m \yi, Vn)(N\Q yi ■ ■ ■ Q yn \I N ) = RHS. (2.30) 

1<2/1<— <Dn<L 

Here Cjv, n is the constant appearing in (I2.27p . In the third equality we used (I2.12p and 
in the forth equality we used the fact that Yli< yi < < Vn <L \Vu • • • > Vn) (z/ij • • • > Vn\ ac ^s as an 
identity in the subspace with n particles. □ 

This is a generalization of the well known duality for the symmetric simple exclusion 
process j9]. The computation of k point correlation functions of (I2.20p is reduced to the k 
particle problem. 

To study the exponential moments of currents, we need a formula when equal. It 

turns out that the quantity can not be written as a summation of only n-particle transition 
probability but as a sum of k particle ones for all k(< n). The result is 

Proposition 3. 

n 

(Q2)t = E^- 1 ) ■ ■ ■ (r n - k+1 -l) ]T £ G({x} k] t\{y} k] 0)(Q m ■ ■ ■ Q Vk ) . 

k=0 l<x\<...<Xf l <x l<J/i<...<yfe<ij 

(2.31) 

To derive this we need a few lemmas. One first shows 
Lemma 4. 

x—l 

(N\Q 2 X = (r - 1) J2( N \QiQ* + ( N \Q* ( 2 - 32 ) 

Proof of Lemma \4\ First one computes 

LHS = C N (0\(X~) N Ql = C N (0\[(X-) N ,Ql) = C N (0\[(X~) N ,Q X ]Q X . (2.33) 
Using Lemma [TJ this becomes 

T N - 1 

cw — r^KQ.x- - q^x^xx-)^-^. (2.34) 

T — 1 
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Rewriting (X ) N 1 Q X = [(X ) N 1 ,Q x ] + Qx(X ) N 1 , applying Lemma [T] again and using 
( ESJ) , it is 



x-1 



t- l /*— - 1 ^-^(o| ( (r - 1) J2 44 ) ( X ~) N ~ 2 + (v\4Q*(x- 



\N-1 



r - I I t - 1 



= (r-^C^^O'^KX^^ + C^^KX-)^ 1 . (2.35) 

3=1 

By (I272TD this is the RHS of (1232]) . □ 

Let us define for a fixed x 

B n = (N\ Qn--'Q*n (2-36) 

l<Xl<-<X n <X 

Then we have 
Lemma 5. 

x _ 1 —n+l i _n 

-B n > y Qj = — ; -B n+ i + £? n (2.37) 

z — ' 1 — r 1 — r 

i=i 

Proof of Lemma First we see 

x n 

Bn Qi = ( U + l ) B n+l + Yl B ^ ( 2 - 38 ) 
3=1 i=l 

where 

B n ,= Yl (N\Qx 1 ---Ql l ---Qx n , (2.39) 

l<a;i<-'-<a;j<3; 

This can be rewritten in terms of B n , B n+ \ only as 

B nii = (t* -VB^ + r^Bn. (2.40) 

This is seen as follows. Using (12.321) . one has 

B nA = (r - l)B n+1 + B n , (2.41) 
B n4 = (r- l)(B nA + ■■■ + B n ^) + i(r - l)5 n+1 + S n . (2.42) 

Suppose f)2.40p is correct for 1,2, — 1. Then one gets (12.401) for i by mathematical 
induction. Plugging (jZHD into (127381) . we get (12737]) . □ 

Now using lemmas I4"f5| it is not difficult to show proposition [3] by mathematical induc- 
tion. 
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3 Step Markov initial condition for the ASEP on Z 



To consider the ASEP on Z, we first put the ASEP with reflective boundaries of size 2L + 1 
on {— L, —L + 1, . . . , L — 1, L}. Then by taking the L — >■ oo limit in ( I2.3ip . we have, for 
the ASEP on Z, 



(Qx)t = ^2(^ n -!)••• (r n - k+1 - 1) 

£ £ G({*}*;t|{y}fc;0)(g sl ---g Mk > . (3.1) 



fc=0 

X 

— oo<xi<...<a;j.<x — oo<j/i <...<j/fc<oo 



Here G({x} k ; t\{y}k', 0) is the transition probability of the ASEP on Z and is defined 
by (J2I9D,(|22DD with 

X 

N x = n r ( 3 - 2 ) 

j=-oo 

Let us assume q > p from now on. For the summation in (13. ip to converge, we assume 
in the sequel of the paper that there are no particles far to the left and there are enough 
many particles far to the right. With this in mind we state the formula as 

Proposition 6. The nth moment of Q x is written as 

n 

(Qx)t = X}( T ™ _ !) " " (r n ~ k+1 - l)c k (3.3) 

k=0 

with 

c k = Yl E G({x} k ;t\{y} k ;0)(Q yi ---Q yk ) . (3.4) 

— oo<xi<...<o;)i ; <a; — oo<yi<...<j/^<oo 

Let N t (x) be the integrated current at bond between sites x and x+1, i.e., the number of 
particles which hop from x + 1 to x minus the number of particles which hop from x to x + 1 
up to time t. Notice that when we consider the initial condition such that rj x = 0, x < 0, 
then N t (x) is the number of particles on sites < x and hence (Q™)t — (r nNt ^ x ') for x < 0. 
Hence the quantity in Proposition |6] is the same as the exponential moment of the current 
of the ASEP. For x > and more general initial conditions, one has to modify the relation 
between Q x and N t (x) to incorporate the initial configuration of particles. 

For the ASEP on Z, the transition probability is written as [T9|l25] 



G({x} k ;t\{y} k ;0)=J2 [ ~'l #i • • ' ^ k A a J] Q^^e^A^ (3.5) 
where S k is a set of all permutations of order k, e(£) = p/£ + qt; — 1 and 

At — sgna — - ■ — — . {6.0) 
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Cr is a contour enclosing the origin anticlockwise with a radius large enough that all the 
poles in A a are included in Cr. In [25J, the contour was taken to be a small one, but one 
can simply use the transformation £ — > l/£ to switch to a large contour. 

Suppose we substitute this representation into (13. 3p . Taking _R > 1, one has > 
1, 1 < i < A; so that the summation over a; can be performed as 

. < 4<T(1) ' ' ' 4(j(fc) " (e.d) -!)••• (e.(D ■ • • - 1) ' 1 j 

By using (13 .5 p . (13 -7p and a formula given in [25] . 

4 8gaff (^(i)-i)-(U)-^)-i) " q IL(i-6) ' ( } 

one arrives at 

Theorem 7. TTie n£/i moment of Q x is written as ([Oj) 



x £ %:i* )0 - (3-9) 

-oo<yi<...<y fc <oo 1,1 ^ k 

In this expression, the dependence on the initial condition is clearly separated. It is straight- 
forward to check whether the summation over y can be taken for a given initial condition. 

From (13.91) one guesses that the distribution function of the current is given by (again 
for a special case where T] x = 0, x < initially) 

F[N t {x) >m] = (-l) w T i ro(m -^ Vt^I * ~ M c fc , x < (3.10) 

fc >o \k-mj r 

where i^) T is r-binomial coefficient difined as {^) T = ^~?i~^~7^^~ ^ I n f ac t following 
the argument in [30J, one sees that ( I3.10p leads to the expression, 

oo fe— 1 

( e AiV t (z)} = ^ r -|fc(fc-l) e Afc "Q^ _ e -A r i) . Cfc _ ( 311 ) 

fc=o j=o 

When A = ralogr, n 6 {1,2,.. .}, the summation over fc is terminated at k = n and this 
reduces to ( 13. 3p . Of course to go in the opposite way from ( 13 .3p to (13. lip , there is a question 
of analytic continuation, an infamous problem in replica theory. We do not discuss it here 
but just mention that it looks natural to expect (13. lip from (13. 3p . 

To further proceed, we need to take the summation in (13. 9p . As explained in [23], there 
are not many examples for which this has been done. Here we will give a generalization 
by observing 
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Lemma 8. Suppose {Q yi ■ ■ -Q yk )o has the form, 



(Qyi " " 'Qy k )o — a fclyi>l Jffl'iCz/i 



Ui-i 



(3.12) 



i=l 



with the convection yo = 0. Here a k does not depend on yi and the function gi,l < i < k 

o |irr converges for |£| large enougi 

{Q yi ---Q yk )o „ TJST^ 9i{Vi 



is assumed to be such that Yly^Lo f^rr converges for |£| /arge enough. Then 



E 

l<yi<...<2/ fe <cx3 



fc oo 

" n yi 

k i=l ?/ 4 =0 



(Ci • • • S^ 



(3.13) 



Proof of Lemma\^ Due to in ( I3.12p . the summation in (I3.9P can be replaced by 

£)i<yi<tti<-<i&- B y shifting Vi^Vi~ 1, the LHS reads 



oo oo 

a <=E e 



E 



9i{Vi) 02(2/2 — 2/1 — 1) fc (y fc - y k _ x - 1) 



SI S2 Sfc 



J/2 + 1 

yi=0y 2 =yi+l 2/fe=2/fe-i+l S1 s2 



(3.14) 



.4 



One further makes a change of variable, yi — >• y± + (y 2 + l) H (?/« + 1), to get the RHS. □ 

Let us consider the initial condition in which there is a particle at the origin and the 
particle occupation on Z + = {1,2,...} is described as a Markov process with the 2x2 
transition matrix, 

(3.15) 

where < p, p: < 1. This means that if the site x(> 1) is empty then the site x + 1 is 
empty (resp. occupied) with probability 1 — p (resp. p) and if the site x(> 1) is occupied 
then the site x + 1 is empty (resp. occupied) with probability 1 — p, (resp. p). We call 
this the step Markov initial condition. When p = pi, each site x > 1 is independent and 
this becomes the step Bernoulli initial conditions [29] except that there is a particle at the 
origin. When the site x is occupied, let us denote by piijn) the the probability that there 
is a particle at site x + L + 1 and that there are m particles on sites from x + 1 up to x + L. 
We set 



(C) = X>HC m 



(3.16) 



m=0 



Notice, if one defines 



Z L {Q = {1\[A 



1 

o C 



|1>= '(0,1), 



(3.17) 



this is written as wl(C) — Z l(C) I Z l(X) ■ 

Let yo = 0, 1 < j/i < 7/2 < • • • < Vk- Suppose there are particles on yi, 1 < i < k and 
that there are rrii particles on sites between + 1 and yi — 1, 1 < % < k. This happens 
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with probability 

k 



Y[Pyi-yi-i-i( m i) ( 3 - 18 ) 

i=i 

due to the Markov properties of the measure. We also have 

k 

Qyi--- Qy k = II T mi+ - +m+i = T ^ k{k+1)+ ^ ik - i+1)mi . (3.19) 



Hence 



i=l m,i=0 
k 

= r i*(*+i)TT„, J^-i+i 



^II^^-H^ 1 ). (3.20) 



1=1 



This is exactly the form in Lemma [8] with gi(y) = w y (r and hence 



E {Qyi ■ ■ -Qy k )o _ lfc(fc+l) TT /o 9 i\ 
" iiz^ te ... a w + i- (3.21) 

Combining Theorem [7] and Lemma [8], we summarize our result as 

Theorem 9. For the step Markov initial condition described around i3.15\) . the nth mo- 
ment of Q x is given by $3.3]) with 



c k = ( _i)* ff i*(»-i) T i*(*-i) r ... r d6 . . . de* n A/ 8 , n f < 



^w+i 

t=l yi=0 



It may still be in general nontrivial to take the summation over yi in f!3.22j) explicitly but 
the point here is that the summations over y^'s are now separated and hence one may be 
able to do asymptotics using this expression with (13.161) . (13. 171) . 

In the parallel way, one can also treat the initial condition in which there is no particle 
at the origin and the particle occupation on Z + = {1,2,...} is described as a Markov 
process with the transition matrix (13.151) . The only difference is that the first product 
p 2/1 _i(mi) is replaced by in (13.181) where p£ > {m) is the the probability that, 

when the site x is empty, there is a particle at site x + L + 1 and that there are m particles 
on sites from x + 1 up to x + L. We also define w^\Q accordingly. Then the net change 
in (I3.22p is that one replaces r5 fe ( fc+1 ) by T^ k ( k ~^ and w yi (r k ) by w y \ > (T h ) respectively. 
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When p = p, the measure (assuming there is no particle at the origin) becomes the 
step Bernoulli initial conditions. For this case, one has 

p L (m)=p( L ) P m (l-p) L - m , (3.23) 

L 

MO = E PL{m)C = P(l - p(l - C)) L - (3-24) 



m=0 



Applying Lemma [8J with = p ,gi(y) = (1 — p(l — C)) V : one finds 

E9i(lli) 1 
m=0 = ^■■•^-l + p(l-r fc -+ 1 )- (3 ' 25) 

The final result is given by 

y- (Qyi---Qy k )o = ik(k-i) fc A 1 (3 26) 

This agrees with the expression given in [27J. 

Next we consider the m-periodic initial conditions in which particles start from the 
sites l,m + 1, 2m + 1, . . . [8]. In this case one has 

Q y = l y>1 r V -^ (3.27) 
where = means the modulo m. Hence 

if. -=C < oo VI- 1 oo TO- 1 oo Vk- 1 

tvi . . tvk 2-^ fV 1 2-^ fV 2 "' 2-^ fVk ' yo.zo) 

-oo<y 1 <...<y k <oo ^ 1 ^ k ai =l ^ s/ 2 =!/l+™» 2 Vk=Vk-i+ m 

yi = l 2/2=1 yfc =l 

This is not exactly the form treated by Lemma |HJ But if we make a change of variable, 
Hi = myi — m + 1, 1 < i < k, then 

%'''^ )0 = T^-Vl^t JJ r (fc-*+i)(fc-*-i-D. (3.29) 
Now one can apply Lemma [HJ with ^(y) = r^ fe+l_1 ^ to obtain 

E (Qyi • • • Qy k )o _ lk(k-l)(t (: yn-lTJ \ (llVh 

-oo<j/i<...<y fc <oo ^ i=l ' 

This agrees with the expression given in [81130]. 
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